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We investigate the two-photon decay rate of a highly excited atomic state which can decay to 
bound states of lower energy via cascade processes. We show that a naive treatment of the process, 
based on the introduction of phenomenological decay rates for the intermediate, resonant states, 
leads to lower-order terms which need to be subtracted in order to obtain the coherent two-photon 
correction to the decay rate. The sum of the lower-order terms is exactly equal to the one-photon 
decay rate of the initial state, provided the naive two-photon decay rates are summed over all 
available two-photon channels. A quantum electrodynamics (QED) treatment of the problem leads 
to an "automatic" subtraction of the lower-order terms. 

PACS numbers: 12.20.-m,12.20.Ds,31.30.jc 



I. INTRODUCTION 

The problem of the calculation of the two-photon decay 
rate from the hydrogen 2S to the 15 state was solved in 
1931 by Maria Goppert-Mayer [1]. Essentially, Goppert- 
Mayer interpreted the two-photon decay in second-order 
time-dependent perturbation theory as a combined tran- 
sition, from 25" to a virtual P state and then, in a sec- 
ond step, to the ground state. However, if one gener- 
alizes the problem trivially, namely to the SS* — > IS* 
two-photon decay, then one inevitably encounters cal- 
culational problems due to the presence of a resonant, 
intermediate, virtual 2P state which leads to a quadratic 
singularity along the photon energy integration contour 
(see Refs. [H US E H f7J ) . Formally, thus, in the limit of 
vanishing decay width of the intermediate resonant 2P 
state, the total two-photon decay rate for 35' — > 15 would 
be infinitely large, provided we apply the formalism em- 
ployed by Goppert-Mayer to the 35 15 decay without 
any modifications or regularizations. 

Of course, this problem has been realized, and three 
possible considerations have been used in order to cir- 
cumvent it: (i) We may observe that the 35 state can 
decay via one-photon, electric-dipole decay (35 2P, 
see p. 266 of Ref. [5]), and we can therefore safely assume 
that the decay rate of 35 will be given by the one-photon 
decay rate 35 2P to an excellent approximation, and 
thus we may entirely neglect the two-photon contribu- 
tion, (ii) Realizing that the virtual 2P state causes prob- 
lems, one can speculate about a possible exclusion of this 
state from the sum over the virtual, intermediate states 
in the propagators [H EHH [3 [S] ■ However, it has been ar- 
gued in Ref. [S] that this procedure is not gauge invariant 
and therefore cannot lead to a consistent solution of the 
problem, (iii) Recently [TOj [11] , it has been pointed out 
that the divergences associated with the quadratic singu- 
larities do not occur if one interprets the two-photon de- 
cay rate as the imaginary part of the two-loop self-energy. 
This observation is in full analogy to the one-loop self- 
energy whose imaginary part gives the one-photon decay 



rate [12]. In the treatment used in Refs. [9j[T0j[TT], one 

obtains expressions for the two-photon decay rate which 
are finite as the regulators (infinitesimal ie terms in the 
propagator denominators) approach zero at the end of 
the calculation, even if we investigate the problematic de- 
cay modes 35 15, 45 15 etc. Besides the analogy 
with the one-photon decay rate/one-photon self-energy 
relations, the treatment used in Refs. [ni[ini[II] has been 
made plausible by field-theoretical arguments and analo- 
gies to radiative corrections to cascade processes in pho- 
ton emission in crossed electric-magnetic fields [13. and 
by analogies with the treatment of quadratic singularities 
in Lamb-shift calculations for highly excited states [T4| . 

However, none of these treatments answer the simple 
question: What do we do about the expressions used by 
Maria Goppert-Mayer, when generalized to the 35 15 
or 45 — > 15 decays in an obvious way? What if we regu- 
larize the divergences due to virtual low- lying P states by 
the most physical, most obvious regularization available, 
namely the total, physical decay rate of those intermedi- 
ate, resonant, P states that we insert into the propaga- 
tor denominators? (This approach has been discussed at 
various places in the literature, e.g., p. 2447 of Ref. [5] 
or p. 4 of Ref. [7].) How can we interpret the result- 
ing expressions from this perhaps nai've, but physically 
intuitive and plausible approach? 

We attempt to answer these questions here. Essen- 
tially, we find that the nai've generalization of the ex- 
pressions used by Goppert-Mayer to the 35 ^ 15 and 
45 — *■ 15 decays is perfectly reasonable, provided we sub- 
tract a lower-order term which is hidden in the expres- 
sions for the 35 — > 15 and 45 15 decays, but absent in 
the 25 15 decays. After the subtraction, the result ob- 
tained using the nai've formalism is in perfect agreement 
with the result obtained for the two-photon decay rate 
from the imaginary part of the two-photon self-energy. 
We reemphasize that the result for the 25 ^ 15 decay 
rate obtained by Goppert-Mayer does not require any 
additional subtractions or modifications. However, one 
should realize that in order to extract physically valid in- 
formation from the naive expressions for the two-photon 
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decay rates from highly excited states, Hke 3S ^ IS 
and/or 45' — > IS* etc., one first has to subtract a lower- 
order term, whose physical interpretation is also dis- 
cussed in the current paper. 

The general outline of this paper is as follows. In 
Sec. |TT] we qualitatively explain the essence of the con- 
siderations needed to fully understand the removal of the 
lower-order terms which are present in the naive expres- 
sion for the two-photon decay rate and thereby anticipate 
a few considerations to be explained in greater detail in 
the following. In Sec. |III[ we then calculate the lower- 
order terms explicitly for the case of the 43 nS decays 
(n = 1, 2, 3), and explain their removal. Conclusions are 
reserved for Sec. HVl 



II. GENERAL DISCUSSION 

We proceed by way of example. Let us remember that, 
e.g., the naive expression for the two-photon decay 45* 
nS reads (see, e.g., p. 2447 of Ref. or p. 4 of Ref. [f ) 



4ck f 



E 



{nS\\x\\uP) {iyP\\x\\AS) {nS\\x\\vP) (vP\\x\\AS) 



E. 



iS 



Eyp — LU 



EnS — E^p 



(1) 



with n 1, 2, 3. We use natural units (?i = c = eg = 1). 
For the reduced matrix elements, we use the conventions 
of Refs. [H]. In the sequel, we use the terms "decay 
rate" and "decay width" synonymously for F, adopting 
the general convention that F is the decay width of a de- 
caying state (Gamow vector, Ref. [17]) with time depen- 
dence exp(— Fi) for the persistence probability. The ex- 
pression Ti^f denotes the two-photon decay rate (naive 

expression) from state \i) to state |/), whereas F-^^ de- 
notes the total one-photon decay rate of state \i) to 
all possible final states (we here use the dipole/long- 
wavelength approximation exclusively in all calculations, 
and we restrict ourselves to a nonrelativistic hydrogenlike 
ion in all concrete calculations) . The QED expression for 
the two-photon decay rate |5J [101 E] reads 



.(2) 



— lim 



4s^„s 277rm2 



J duJUJ^ {E4S - Eni 



{nS\\x\\iyP) {iyP\\x\\AS) 



E. 
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E,, 



nS\\x\\i^P) {iyP\\x\\AS) 

EnS - E„p + UJ + lt 



(2) 



(2) 

The QED two-photon decay rate is denoted F^j^^ in the 
current work in contrast to the naive expression, which 
carries no superscript. 

In order to understand the difference of the naive and 
the QED expression, we should first observe that the de- 
cay width (even the one-photon width) in the denom- 
inator of the naive decay width is parametrically sup- 
pressed with respect to the energy differences of atomic 
states [f[,p a{Za)'^ vs. £^45 — E^p ^ (Za)^ in units of 
the electron mass], and we can therefore apply the (Za)- 
expansion to the evaluation of the expression ([T|) . Specif- 
ically, we evaluate ([ij in the limit T^jp £'45 — E^p, 

which is justified because the decay width F^^p is of higher 
order in the (ZQ;)-expansion than the energy difference 
E4S — E,jp. Physically speaking, the approximation is 
justified because the lifetime of a typical atomic level is 
long on the time scale of the oscillation period of radia- 
tion emitted in a typical atomic transition. We can thus 
expand the naive expression for the two-photon decay 
width of highly excited levels for small decay widths of 
the virtual states, keeping the physical value of the decay 
width at all stages of the calculation. 

The first observation to be made is that physically, 
in the limit of a small decay width of the intermediate 
state f[,p <C Eis — E^p, we can intuitively assume that 
the entire decay rate of the initial state will be due to 
one-photon decay from the initial state to the intermedi- 
ate, resonant states (i.e., to exactly those virtual states 
with intermediate energy between the initial and the fi- 
nal state of the two-photon process). Indeed, once the 
system has decayed to one of the intermediate state, it 
is "stuck there" in view of the small decay width (long 
lifetime) of that intermediate state. Intuitively, we would 
thus expect that the naive expression for the two-photon 
decay rate might contain, as a lower-order term in the 
(Za)-expansion, the one-photon decay rate of the initial 
state. Below, we show by an explicit calculation that this 
is indeed the case, and that the naive expression for the 
two-photon decay rate simply contains the one-photon 
decay rate as a lower-order term. 

This slightly oversimplified statement is shown to hold 
provided one-photon cascade processes through interme- 
diate, resonant states are possible for the chosen initial 
state under investigation, and provided we sum over all 
open cascade channels. (The first condition is not ful- 
filled, e.g., for the 25* — > 15 decay, where no intermediate, 
resonant states are available and there are no lower-order 
terms to subtract.) 

In terms of the (Za)-expansion, the lower-order term is 
calculated to be of the order of order a{Za)'^ as opposed 
to a^^Za)^ (the latter would be the expected order-of- 
magnitude for the two-photon decay rate). The remain- 
ing term from the naive two-photon decay, i.e. the ex- 
pression left over after the subtraction of the term of or- 
der a{Za)'^, then is the coherent two-photon decay rate, 
and this left-over term is equivalent to the prediction 
from QED theory. 
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A final word on Eq. (ITJ is in order. It has been ques- 
tioned in the literature (e.g., in Ref. [7]), whether one 
should use the total or a differential decay rate for the 
propagator denominators in this case. The most physical 
and most intuitive prescription is to use the total decay 
rate of every single intermediate state in order to regu- 
larize the quadratic singularity. Indeed, we find here that 
a very clear interpretation of the lower-order terms can 
be given provided the total decay width of the interme- 
diate, virtual states is used as a regulator. Because the 
total decay rate can be approximated very well by the 
one-photon decay rate for all the virtual states, we use 
the total one-photon decay rate r[^p as a regulator in the 
denominators of Eq. (jlj . One would have to change this 
prescription only if the virtual 25* state were present as a 
virtual states because it is metastable and decays primar- 
ily via two-photon decay (Fj^ = in the nonrelativistic 
dipole approximation), but the 25' states is not present 
as a virtual state in our calculations. 



rate) corresponds to a natural quantum mechanical ad- 
dition of amplitudes for the two photons to be emitted 
(in whichever sequence) before taking the square. In this 
case, for vanishing regulator e — > (here, e corresponds 
to the regularizing decay width) , we obtain divergent in- 
tegrals of the form 



a — oj + le 



-0 TT 



1 



e a(a — 1) 
u 

for < a < 1. Combining ([3| and we have 



+ 0{e'), (4) 



duj 



a — uj + le 



^ TT 

e 



a — OJ + le 



A natural generalization of this formula, which will be a 
key to our further considerations, reads 



III. CONCRETE CALCULATION 

First, we recall the mathematical mechanism by which 
the quadratic singularities are regularized according to 
Refs. [HI [ini E]. We use Feynman's ie prescription for 
the propagator denominators and consider the following 
model integral [see Eq. (16) of Ref . 



lim Re 



dcij 



1 



a — UJ + ie 



1 



a{a — 1) 



(3) 



Here, e (0, 1) is a scaled variable that corresponds to a 
scaled photon energy (scaled relative to the entire photon 
energy interval for the decay process) . 

The naive regularization method (insertion of a decay 



fia) 



■ le 



duj 



a — a; + ie 



+ (5) 



with /(w) e R. Equation ^ follows from ^ and Q 
and from the observation that the two integrands differ 
only in a small region about w ^ a. Further illustra- 
tive remarks regarding the model examples are given in 
Appendix [A] 

We now investigate the naive expression for the total 
two-photon decay width of the 45* state, which is the sum 
of the 45 -> 15, 45 ^ 25, and 45 ^ 35 channels, 



Eis-E^s 



4a^ 
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E 
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E 



(25||f||^P)(t.P||f||45) (25||f||^F)(z.F||af||45) \ 



^ yEis — E^p — + I i r[,p E2S — E^p + + ^ i t')^p) 



4a" 



Eas—E^s 



277rTO" 



dww^ {Eis - E3,s - Lof 



E 



{iS\\x\\uP) {vP\\x\\AS) {iS\\x\\vP) {uP\\x\\AS) 



T' \E^s — Ejjp — o; + i i r[,p E'is — E^p + cj + ^ i r[,py 



(6) 



From the sum(s) over we can now single out those terms which generate the double poles in the photon energy 
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integration in the limit Fj^p 0, and a remainder term TZ, which contains at most simple poles in that limit, 
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dio {Eis - Eis - ujf 



4a" 
277rm2 
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Eis — E2P — uj + \ i T^2P Eis — E2P + a; + i i 



(15||f||3P)(3PP||45) , {lS\\x\\3P) {3P\\x\\4S) 
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4q" 
27W 
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(25||f||3P)(3PP||45) , (25||x||3P)(3P||f||45) 



i?4S — E3P — UJ + Fgp i?2S — E^p + a; + ^ i Fgp 



1 ;-r(i) 



n. (7) 



Note that the last term in Eq. (|6|, which corresponds to 
the 45* 3S decay, does not contain any terms which 
could possibly lead to cascade processes, and is thus en- 
tirely contained in TZ. We now use Eq. ^ with the iden- 
tification e = I F^p with n = 2, 3 as appropriate for 
the treatment of the double poles. For the simple poles, 
including those contained in the remainder term TZ, we 
only need a principal-value regularization according to 
the model integral 

Re f\^-[i^=P.Y. [\,Ii^+o{e). (8) 
Jo a-uj + i€ Jq a-uj 

again with the identification e = ^F^p, as appropri- 
ate. The decay width regulators in the denominators 
of Eq. m provide for the imaginary parts. We can ig- 



nore terms of order e in Eq. (|8]) which would otherwise 
lead to corrections to the decay rate of order a^{Za)^. 
Our self-explanatory notation for the partial one-photon 
decay rates is Ti^j for the channel from level \i) to \f). 
Using the formulas 

^llLrnP = ^ " (Za)' (P„s - E,,pf {nS\\x\\mPf , 

^nLrnS = ^ " (^«)' " E^pf {nS \ \x\ \mP)\ 

(9) 

together with ^ and (|8|, we can finally reformulate the 
nai've two-photon decay rate F45 as 
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2 2P 
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2 3P 


(2511 
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27^"^' J \E4s - E:iP - UJ + \iT^^)> E2S ^ E^p + uj + llT^^p ^ 
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(10) 
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We can now drastically simplify those terms in the above expression which are free from photon energy integrals, 

(1) ^ 2P— IS-L 4S~*2P (1) ^ 3P— 25-^ 4S^3P (i) ^ 3P— IS"^ 4S^3P 

^ 2P ^ 3P ^ 3P 



-(1) _L _ 

4S^2P + (1) 
^ 3P 



3P^25 



-(1) 



3P^1S i ^ 4S^3P 



-(1) 



4S— 2P ^ ^ 4S^3P 



= r 



(1) 
4S ' 



(11) 



where we have used — r2P^i5- This consideration implies that the term n/e in Eq. (jlj, when identified as 
e = ^ r^^p (with n = 2,3) and applied in the model example in Eq. (|5|, generates the full one-photon width F^^, 
provided we sum over all open two-photon channels. Thus, we can write for the total nai've expression of the decay 
rate, which is composed of the sum of the 45" IS, iS 2S, and 45" 3S decays. 
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(12) 



This equation simply means that the naively regularized 

two-photon decay rate can be written as the sum 

of three terms: (i) the total one-photon decay rate F^g 

(which consitutes the spurious lower-order term), (ii) the 

(2) 

QED two-photon correction F^^ and (iii) higher-order 
terms 0{a'^{Za)^). Note that the second, third and forth 
term in the expression following the first equal sign of 

^ 5 I'l/P' ^^'^ J^^t the 



Eq. (12 1, with the identification e = ;^FJ 



QED expressions for the two-photon decay widths corre- 
sponding to 45" IS, AS 2S, and 45* — > 3S, respec- 
tively. As shown in Refs. [9j[10j[TT], the two-photon decay 

rates F^^^^^^ are of order a^(Za)^, as they should be. 
We also observe that the QED treatment [5J [101 Ej thus 
eliminates the lower-order terms right from the start. 

The generalization to, e.g., an initial 55" state is 
straightforward, but one has to take into account the 
55 — + 4P — > 31? cascades. We will refrain from dis- 
cussing the generalization further, as the principle of the 
calculation should be sufficiently addressed by the exam- 
ple case of the 45* state. 



IV. CONCLUSION 



The central result of this paper is contained in Eq. ( 12 1 



For the initial 45* state, it is shown that the coherent two- 
photon contribution to the decay rate, as given by a QED 



treatment, is obtained from the naive expression for the 
two-photon decay width, given in Eq. ([ij, after the sub- 
traction of lower-order terms whose sum over all open 
channels is exactly equal to the total one-photon decay 
width of the initial state. Note, in particular, that the 
QED treatment eliminates the lower-order terms right 
from the start [SJ [101 HI]- This result could be antic- 
ipated (see Sec. |ll]) based on a physical consideration: 
The (ZQ;)-expansion of the decay rate implies an expan- 
sion for small decay widths of the intermediate, resonant 
states; this expansion, in turn, physically implies a long 
lifetime for the intermediate, resonant states (on the time 
scale given by the oscillation period of radiation emit- 
ted in a typical atomic transition). This means that in 
a first approximation, the virtual resonant states "trap" 
the electrons after a single photon has been emitted. This 
again implies that the lower-order term contained in the 
naive expression for the two-photon decay rate should be 
equal to the one-photon decay rate of the initial state, 
provided we sum over all open channels, and provided all 
possible one-photon decays of the initial state can par- 
ticipate in such two-photon cascades. The term left over 
after the subtraction of the lower-order terms is the two- 
photon correction to the decay rate, and it is of order 
a^{Za)^ in units of the electron mass, as it should be. 

These considerations allow us to give a very clear phys- 
ical interpretation of the lower-order terms which are 
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to be expected in the treatment of a two-photon decay 
process with virtual resonant states: Because the addi- 
tion of the probabihty amphtudes in the naive treatment 
docs not lead to a clear separation of the sequential two- 
photon emission via the intermediate, resonant states 
from the coherent two-photon contribution to the decay 
rate, lower-order terms have to be expected, and in view 
of the above arguments, these have to be exactly equal 
to the one-photon decay rate of the initial state. Let us 
recall what a treatment of the total decay rate would im- 
ply for the decay width of, say, the AS state, if we were to 
add the naive expression for the two-photon decay rate 
to the well-established result for the one-photon decay 
rate [8] of that state. In this case, because of the presence 
of the lower-order term in the two-photon decay width, 
we would double-count the one-photon decay rate of the 
45' state and obtain a result for the total decay width 
which would be twice as large as the established value. 
This procedure cannot be consistent, and a subtraction 
is definitely required, as outlined here. 

Note that spurious lower-order terms are a recurrent 
theme in a number of QED calculations. E.g., in the 
first, non-relativistic treatment of the bound-electron 
self-energy [18^, the result was found to diverge only log- 
arithmically with the energy of the virtual photon, but 
only after a more severe linear divergence was identified 
as a low-energy part of the mass renormalization asso- 
ciated with the low-photon-energy contribution to the 
Lamb shift and subtracted (see also Refs. [TH UHl 1^). 
Similarly, in Ref. [22 the authors had to subtract 
several spurious lower-order terms in order to obtain 
the a{Za)^ correction to the Lamb shift from the two- 
Coulomb-vertex forward scattering amplitude [see espe- 
cially Eqs. (36) and (53) of Rcf . [23 ]. Eventually, these 
spurious terms could be shown to cancel against com- 
pensating terms from the one-Coulomb-vertex term, but 
if the authors had considered only the two-vertex con- 
tribution to the Lamb shift, then they would have had 
to carry out a subtraction, and this is indeed the spirit 
in which most modern Lamb shift calculations are done 
[see also the remarks after Eq. (32) of Ref. [23], where 
several lower-order subtraction terms are identified and 
subtracted, albeit in a different context]. While the result 
reported here could have been found immediately after 
the initial treatment due to Goppert-Mayer [T], concepts 
developed later in the context of the theory of renormal- 
ization have inspired us to look for possible subtraction 
terms. Note that the two-photon decay rate can natu- 
rally be interpreted as the imaginary part of the two-loop 
self-energy [TOj and therefore also constitutes an "energy 
shift," albeit an "imaginary one," or an "imaginary part 
thereof." The additional subtraction term found here can 
thus be interpreted as an additional renormalization of 
the imaginary part of the bound-state energy which be- 
comes necessary if the naive method of regularization is 
used for the resonant intermediate states. 

Summarizing the results of our investigations, we con- 
clude that physically sensible lower-order terms are ob- 



tained from the naive expressions for the two-photon de- 
cay rates when the propagator denominators of the inter- 
mediate, resonant states are regularized by assigning to 
them the total decay widths of those virtual state. This 
is in agreement with physical intuition and reassuring. A 
final remark on gauge invariance: For an infinitesimal ie 
in the propagator denominators |10l 111] , all expressions 
are gauge invariant with respect to length and velocity 
gauges, as shown in Ref. [S]. For finite, small decay rates 
acting as regulators, gauge invariance, strictly speak- 
ing, does not hold. However, the gauge-noninvariance 
is shifted to higher-order terms ^ a^{Za)^, which are 
not relevant from a phenomenological point of view, and 
thus, with the correct prescription for the subtraction of 
lower-order terms, the result for the two-photon correc- 
tion to the decay rate can indeed be obtained from the 
naive expression regularized using the decay widths of 
the virtual, resonant, intermediate states. 
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APPENDIX A: FURTHER EXPLANATORY 
REMARKS 



We would like to illustrate the calculation of the two- 
photon decay rate from highly excited states by some 
further remarks. 

At the heart of the calculation lies the model example 
in ^ which gives the basic regularization of the two- 
photon decay rate by the ie prescription. Surprisingly, 
the regulation can also be made plausible, by a principal- 
value integration, although this is not obvious in the case 
of a quadratic singularity. In the related treatment of cas- 
cade processes in field theory [see Ref. [13] and Eq. (6.20) 
ff. of Ref. [23]] and of the corresponding double poles, the 
principal-value integration leads to the result [see also 
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Eq. (32) of Ref. 0] 



duj P.V. 



= lim / dcj P.V. I ( P.V 

Jo V a — uj + T] J \ a - uj 



lim - / du} P.V. P.V. 

»?^o T] Jq \ a — UJ a — Lo + rj 



lim — 
»?^0 ?; 



In 



1 - a 



In 



1 — a — 77 



a{a — 1) 



(Al) 



The results on the right-hand sides of Eqs. (|3| and ( Al ) 
are identifcal. The term n/e accounts for the difference 
among the two results given in Eqs. (|3| and Q; here 
this term is identified as a spurious lower-order term as- 
sociated with one-photon decay to the intermediate reso- 
nant state(s). After its subtraction, the two results using 
either the QED prescription ([3| or the naive prescrip- 
tion Q are in agreement. 



1 


2 

-f 


1 


2 


2tt 


loq — Lu + ie 




LOi — UJ + ie 




e 



The subtraction procedure can also be illustrated 
graphically, as in Fig. [l] by plotting a "model function" 
/(w,^) in terms of its two arguments, which we assume 
to be scaled dimensionless parameters G (Oil)- For 
^ = 1, we plot the naively regularized model function 



f{uJ,^^l)^ 

UJQ — UJ -\~ iti LUl — -r It 

(A2) 

illustrating the naive regularization. This function has 
two maxima near ujq — 0.3 and uji — 1 — ujq which ap- 
proximate the peaks in the two-photon spectrum under 
the presence of a virtual resonant state which is displaced 
from the initial state of the two-photon process by the 
energy ujq — 0.3. We use a numerically small regulator 
value e = 0.01 in accordance with the basic assumptions 
on which our derivation is based; if e were large, then, 
e.g., argument regarding the "trapping" of the electron 
in the intermediate states as given in Sec. |IV| would break 
down. The term — 27r/e is subtracted in order to ensure 
that the function f{uj,^ — 1) has the same integral over 
u! G (0,1) as the QED regularized expression 




/(c^,C = 0)=Re 



1 



1 



{ujQ — UJ + ie)^ {uji — UJ + ie)^ j ' 

(A3) 

which has two minima near ujq and uji and illustrates the 
ie prescription. The interpolation 



(O 



(A4) 



FIG. 1: (Color online) Graphical illustration of the subtrac- 
tion procedure. The explanation is in the text. 



illustrates that in view of the negative tilt of f{uj, ^) along 

the 5 axis, the integral duj f {uj , S,) , i.e., the "decay 
rate," is independent of ^ and thus the same in both 
regularizations (^ = versus ^ = 1). 
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